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ABSTRACT 

Formation of droplets due to breakup of liquid column is important research. In the present study break up of 
planer liquid sheet has been investigated using linear stability analysis. In this study both sinuous and varicose modes of 
disturbance are used. As reported by the earlier studies linear analysis of sinuous mode will not create any breakup of the 
liquid sheet but linear analysis of varicose mode of disturbance can create breakup of the liquid sheet. Also breakup of the 
liquid sheet is very much nonlinear especially at the time of break up. So a nonlinear analysis is required for the complete 
understanding of the underlying physics of the breakup process. But which will be the dominant mode of break up process 
and what will be the frequency and growth rate of fundamental mode of perturbation can be assessed by the linear analysis. 
So in this study temporal linear stability analysis has been done for both sinuous and varicose modes of disturbance to 
calculate the growth rate and frequency. The result shows very good agreement with the earlier studies. Parametric studies 
to understand the influence of fluid velocity ratio, density ratio and surface tension on the growth rate have also been 
investigated. 
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Table 1 



h 


Half thickness of liquid sheet (m) 


k 


Wave number (dimensionless) 


K 


Cut-off wave number (dimensionless) 


"■max 


Critical wave number, corresponding to maximum growth rate (dimensionless) 


L 


Non-dimensional break-up length 


I 


Break-up length of the liquid sheet 


U 


Gas-to-liquid velocity ratio, m/m, 


u, 


Velocity of the liquid sheet 


u g 


Velocity of the gaseous phase 


We 


Weber number (= Pjiifh 1 a ) (dimensionless) 


X 


Axial coordinate (m) 


y 


Transverse coordinate (m) 




Greek Symbols 


a 


Angular frequency (dimensionless) 


P 


Growth rate (dimensionless) 


ftmax 


Maximum growth rate (dimensionless) 


<P 


Dimensionless velocity potential 


T](x,t) 


Instantaneous and local transverse coordinate of interface (dimensionless) 


no 


Initial disturbance amplitude (dimensionless) 


¥ 


Function representing temporal variation of r\ 


CO 


Temporal frequency 


Pi 


Liquid density 
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Table 1 Contd., 



n 
F 


frn^-to-limiirl Hpnsiitv ratio 


P g 


Gas density (kg/m 3 ) 


a 


Surface tension coefficient 


U 


Viscosity of gas 


M, 


Viscosity of liquid 


8 


Gas phase 


1 


Top interface (j = 1 or 2) 


2 


Bottom interface 


I 


Liquid phase 


n 


nth order 


,t 


dldt 


,x 


dldx 


,xx 


8 2 ldx 2 


,xy 


8 2 1 Bxdy 


>y 


dldy 


>yy 


d 2 ldy 2 



INTRODUCTION 

Liquid sprays are encountered in nature and in a wide range of science and engineering applications. It generally 
consists of different sizes of droplets dispersed in a gaseous medium. In nature, droplets can be found as rain, dew, fog, 
clouds and waterfall mist. Mechanical devices like nozzles or atomizers also produce spray. In many industrial applications 
liquid is used as a spray instead of bulk liquid such as spray drying, spray cooling, spray combustion, spray deposition, 
thermal spray, surface treatment, spray inhalation, crop spray, spray painting etc. So spray has become a subject of interest 
in many engineering disciplines like mechanical, aerospace, chemical, metallurgy, materials, pharmaceutical, food 
processing, and agriculture, meteorology and powder generation. Atomization of the bulk liquid can be achieved in 
different ways, such as aerodynamically, mechanically, ultrasonically, or electro statically. Depending upon the 
atomization process droplet size in the spray will also vary. In most of the atomizers, droplet size generally varies from a 
few microns to around 800 - 1000 microns. 

For sprays in agriculture and fire suppression system, the droplets should be sufficiently large to have enough 
momentum so that they could not be carried away by the wind. On the other hand, bronchial sprays require very small 
droplets of the order of few microns, so that the droplets can pass through the small airways and reach the lung. Therefore 
disintegration of liquid sheets and leading to generation of fine droplets has been extensively investigated owing to its 
importance in atomization and spray formation. Excellent reviews on the subject are available in Lin and Reitz 
[1], Mehring and Sirignano [2] and Lin [3]. The injector configuration can be idealized as a liquid sheet sandwiched 
between two gas streams. Annular liquid sheets, which are the configuration of the liquid sheet emerging from most of the 
atomizers, can also be reasonably simplified to a planar liquid sheet sandwiched between two planar gaseous layers if the 
sheet thickness is much smaller than the mean radius of the annulus. 

Earlier studies reported that, though the fundamental sinuous mode generally grows at a much faster rate than the 
varicose mode over a wide range of parameters, still this mode alone cannot cause breakup as the distance between the two 
surfaces are preserved for linear stability analysis. It is the varicose mode, which in spite of low growth rate, causes the 
instability by bringing the surfaces closer. Mitra et al. [4] Nath et al [5] investigated the linear stability of a planar liquid 
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sheet subjected to both sinuous and varicose disturbance. The sinuous mode was observed to be responsible for large 
amplitude deformation while the varicose mode caused the breakup. 

The objective of the present work is to investigate the effect of velocity ratio, density ratio and surface tension on 
growth rate of disturbance for both the sinuous and varicose mode of disturbance. 

I 

Nozzle — 





Ligaments 



1 



/ o 

Droplets 

Figure 1: Schematic Diagram of a Planer Sheet Breakup Leading to Spray Formation Ref [7] 
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Figure 2: Schematic of the Liquid Sheet Interface for (a) Sinuous 
(b) Varicose Mode of Disturbances 
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MATHEMATICAL FORMULATIONS 

Figures 2(a) and 2(b) shows the schematic of a planner liquid sheet surrounding with two gas streams moving 
with non zero velocities and the gas -liquid interfaces for sinuous and varicose mode of disturbances respectively. 
In the undisturbed state, one -dimensional flow is assumed and the velocity of the liquid sheet is given by u, while that of 

each gas medium is given byu g . The densities of gas and liquid layers are given by p and p l respectively. Froude 

number for typical atomizer conditions are generally being high, the effect of gravity has been neglected for the analysis. 
Cousin and Domouchel [6] have shown that for high Weber number (We* /Re, <10 3 ) the effect of viscosity on the 

stability of a planar liquid sheet can also be neglected. 

For the convenience of analysis, all the physical parameters are non-dimensionalized. Length is 
non-dimensionalized by half-sheet thickness h, time by the onvection time h/ul where ul is the velocity of the undisturbed 
liquid sheet, and density by liquid density pl 

When the undisturbed flow, described above, is perturbed by a small disturbance, the positions of the two 
interfaces are described as 

y(x,t)={-iy +l +rjj(x,t) (1) 

In the above equation, j = 1 and 2 denote the upper and the lower interface respectively. The resulting velocity 
field must satisfy the equation for conservation of mass. The interfaces must satisfy the kinematic and the dynamic 
boundary conditions. The surface deformations ^(x,?) must satisfy the above governing equations and boundary 

conditions. The solution to the above is obtained using regular perturbation theory with the initial (small) disturbance 
amplitude j] as the perturbation parameter. The surface deformation is expanded in a power series of r/ as 

. . 00 

r lj {x,t)=Tri%r ljn {x,t) (2) 

n=\ 

Following the analysis of Ibrahim and Jog [9] and Nath et al. [5], the present study will focus on the first order 
solutions only. The governing equations and the boundary conditions for the first order linear stability analysis are as 
follows. 

Governing Equations 

fa\,xx + <t>ll,yy =0, -l + 7 l2 <y<l + r ll (3) 
tgjXxx + tgjXyy = 0; 

-H< y<l + 7] 2 ,-l + 77i ^ y^ H ( 4 ) 
Kinematic Boundary Conditions 

^/i.v-^v-^u =0 (5) 

^%,->7,u-^u = (6) 
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Dynamic Boundary Conditions 

(-iy 

P^u - <f>iu + P L %, x ~ = 77 J.™ (7) 

Solution Procedure 
Solution for Sinuous Mode 

Two major modes of disturbance are generally considered in stability analysis of liquid sheet. Sinuous (when the 
disturbances on the two surfaces of the liquid sheet are in phase) and varicose (when the disturbances on the two surfaces 
have a phase difference of 180°). The two modes of disturbance are shown schematically in Figure 2. 

If the initial amplitude of disturbance is ^° and wave number is k, then the profile of the perturbed surface can be 
described in linear stability analysis by the function 

tjj (x,t) = 7j \y/ jX (t)exp(ikx) + ^7 jl (t)exp(-ikx)j (8) 
Here, y7 n (f) denotes the complex conjugate of y/ jl (f ) . 

Substitution of Eq. (8) in the kinematic boundary conditions, Eqs. (5) and (6) and using the governing equations, 
Eqs. (3) and (4), the velocity potentials for the upper gas, lower gas and liquid sheet are obtained. Substituting the 
expressions of the velocity potentials and the surface deformation in the dynamic boundary condition, Eq.(7), the following 
dispersion equation is obtained: 



l - [tanh (k) + p] y/ ]X {( + 2i [tanh (k) + pU] ys jl{ 

(9) 

M = ° 



9 k 2 

ktanh(k) + pkU 

W We 



¥ iiv) 

The solution of Eq. (9) gives the functional form of 1 

It may be noted that for temporal stability analysis, the wave number k is real while the temporal frequency, co l is 
complex, given by co l =a i + //?, , where a x and /?, represent the angular frequency and the growth rate respectively. 
The initial conditions used to obtain the solution are as follows. 

7] ji (x, 0) = cos(fcc) 

(10) 

r/ji f (x,0) = — OTj sin(£x) 

Applying the above initial conditions, the following solution is obtained: 

y/ ji (/) = ^■cosh(/?jf)exp(;a 1 f) (1 1) 

The growth rate and the angular frequency are given by 

lclp(l-U) 2 tanh(k)- — (p + tanh(fc))l 2 

A = — ^ — (12) 

p + tanh(£J 
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k(pU + Uah(k)) 

and ccj = (13) 

p + tanh(&) 

respectively. 

Solution for Varicose Mode 

The analysis for the varicose mode is similar to that of the sinuous mode. However, as shown in Figure 2(b), the 
deformation of the liquid sheet is symmetric about the mid-plane of the liquid sheet. Thus, only one half of the domain 
needs to be computed using symmetry condition at the central axis (y = 0). The symmetry conditions require that the 
normal velocity and the normal gradient of the axial velocity are zero. In terms of the velocity potential of the liquid, these 
conditions become 

<t>l,y=<t>l,xy=° (W) 

Proceeding similarly as in the case of sinuous mode, the first order dispersion relation, growth rate and angular 
frequency are obtained as 

j [coth (k)-p] y/ jU , + 2i [coth (k) + pU~\ y/ jU 

(15) 

coth(k) + pU 2 



-k 



We 



¥jl =0 



k\p(l-ufcoth(k)- — [p + co±(k))\ 2 

& = T^k\ — (16) 

kipU + coMk)) 

a= \L \1L (17) 

p + coth(£:) 

Fractional powers are preferred to root signs. The solidus ( / ) should be used instead of the horizontal line for 
fractions whenever possible. Consecutive numbers to identity mathematical expressions should be enclosed in parentheses. 
Refer to equations in the text as "Eq.(l),"etc, or "Equation(l)," etc., at the beginning of a sentence. All symbols should be 
in italic letters. 4. 

RESULTS AND DISCUSSIONS 

At all wave numbers and for U < 1, the growth rate for the sinuous mode is at least one order of magnitude higher 
than that for the varicose mode. From Figure 3 it is also clear that for the stability curves for U > 1 and U < 1 the cut-off 
frequency does not depend as strongly on the velocity ratio for U < 1 as in the case of U > 1. Figure 3 also shows the 

variation of maximum growth rate, ^ max with velocity ratio for both sinuous and varicose modes at different Weber 

numbers. The figure shows that for U < 1, ^ max remain invariant with U except for the region close to U ~ 1. 

For values of U close to 1, the growth rate sharply decreases for both sinuous and varicose modes. In this 
formulation, the instability is primarily driven by the velocity difference across the interface. Hence, for velocity ratios 
close to 1, this driving factor gets weakened, leading to small growth rates. From the figure, it is also observed that for all 
Weber numbers and U < 1, the growth rate is higher for sinuous mode. 
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(a) We=150 rho=.001 and U<1, 
Sinuous Mode of Disturbance 
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(c) We=150 rho=.001 and U<1, 
Varicose Mode of Disturbance 
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Varicose Mode of Disturbance 



Figure 3: Variation of Growth Rate with Wave Number for Different Velocity Ratio 
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Figure 4: Variation of Maximum Growth Rate with Velocity Ratio for Different Weber Numbers 

For U > 1, both these parameters increase with velocity ratio and at very high velocity ratios, the sinuous and the 
varicose modes have similar maximum growth rate. The range of velocity ratios for which the two modes show similar 
instability characteristics increases with increasing Weber number. Figure 4 shows the variation maximum growth rate 
with velocity ratio for different Weber number. The growth rate is higher for sinuous mode than for varicose mode, but the 
difference decreases with increase in Weber number. 
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CONCLUSIONS 

A linear model has been developed for predicting the growth rate of breakup for a planar liquid sheet sandwiched 
between two moving gas streams. The velocity of the gas streams is taken to be equal but different from that of the liquid 
sheet. Closed form solutions for the gas and liquid velocity fields and the interface shape are obtained assuming inviscid, 
irrotational flow for both the phases. The effects of gas velocity, both higher and lower than the liquid velocity, are 
investigated for a wide range of Weber number .The major observations of the investigation is for U < 1, the growth rate is 
higher for sinuous mode at all Weber numbers. However, for U > 1, the growth rates for sinuous and varicose modes 
become similar at high velocity ratios. 
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